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Abstract 

There exists a positive function ip{t) on t > 0, with fast decay at 
infinity, such that for every measurable set Q in the Euchdean space 
and R > 0, there exist entire functions A (x) and B (x) of exponen- 
tial type R, satisfying A{x) < xn{x) < B{x) and \B{x) — ^ 
{R dist {x , dQ)) . This leads to Erdos Turan estimates for discrep- 
ancy of point set distributions in the multi dimensional torus. Analo- 
gous results hold for approximations by eigenfunctions of differential 
operators and discrepancy on compact manifolds. 

An infinite sequence of points {xj}^^ is uniformly distributed in the 
interval < a; < 1 if, for every < a < 6 < 1, 



lim <^ m ^y^Xia,b] (xj) \ = b - a. 



A quantitative measure of the irregularity of distribution of the points 
{^j}f=i is given by the discrepancy 



sup 

0<a<fe<l 



{b-a)-m ^^X[a,b] {Xj] 



Upper bounds for the discrepancy are useful because they lead to upper 
bounds for the approximation of integrals by Riemann sums. A well known 
criterion due to H.Weyl states that a sequence is uniformly distributed if and 
only if for every k ^ 0, 

r m ~\ 
lim < m^^y exp {2'Kikxj) ^ = 0. 
K i=i ) 
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An upper bound for discrepancy is given by the classical inequality of 
P.Erdos and P.Turan: 



sup 

0<a<6<l 



(6 - a) - m ^^X[a,b] {xj] 



< 



n=i,2,... 1 n + 1 \A; n + lj 



m ""^^^ exp {2nikxj) 



A proof of the above inequality relies on approximations from above and 
below of the characteristic functions Xla,b] (x) by trigonometric polynomials 
P{x) = X]fc=-n -^(^) ^^P (2vriA;x), so that 



{b-a)-m ^^X[a,b] (Xj) ^ (b 



a] — m 



= (6-a)-P(0)- lm~^Yexp{2mkxj)\P{k). 

l<|fc|<ri V j=l ) 

See [21 [71 [15]. A deep study of the approximation of the characteris- 
tic function of an interval by trigonometric polynomials has been done by 
A.Beurling and A.Selberg, who proved that for every < a < 6 < 1 and 
n = 0,1, 2,..., there exist trigonometric polynomials P± {x) of degree n with 

P~ (x) < X[a,b] (x) < P+ (x) , 

{x) - X[a,b] {x)\dx = l/{n + l). 

See e.g. ^j. Similar extremal problems have been considered in f]3], with 
precise estimates on the approximation, in —oo < x < +oo with measure 
Ixl^*^^^ dx, of the function sgn [x) by functions of finite exponential type. A 
radialization of these functions then yields an analog of Selberg polynomials 
for approximation of characteristic functions of mult i- dimensional balls and 
this has been applied to Erdos Turan estimates of discrepancy. See [SI [131 E], 
and also [5j for a generalization to boxes. 

Here we look for a geometric analog of Erdos and Turan results in a very 
general framework, where the intervals in the torus are replaced by arbitrary 
measurable sets on a manifold, and trigonometric polynomials are replaced by 
finite linear combinations of eigenfunctions of the Laplace Beltrami operator. 
In spite of this generality, the techniques are rather simple, and they may be 
of some interest even in the classical Euclidean setting. Moreover, the results 
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are optimal, up to the constants involved. In fact we shall try to pay special 
attention to these constants, which are quite explicit although not optimal. 
In this sense we acknowledge that we do not match the beauty of Beurling 
and Selberg results, which stems precisely in their extremal properties. 

The plan of the paper is as follows. The first section is devoted to approx- 
imations from above and below of characteristic functions by entire functions 
of exponential type or, in the periodic setting, approximations by trigono- 
metric polynomials. The main result of this section is the following: 

Theorem 1 There exists a positive function il){t) on t > with fast decay 
at infinity, ip{t) < c(a) (1 + for every a, such that for every measurable 
set Q in the Euclidean space M'^ and R > 0, there exist entire functions A{x) 
and B{x) of exponential type R, satisfying 

A{x) < xn{x) < B{x), \B{x) - A{x)\ ^ i; {Rdist {x,dn)) . 

Roughly speaking, the approximation \B{x) — A{x)\ is essentially 1 at 
points with distances 1/R from the boundary of Q, while \B{x) — A{x)\ is 
essentially at larger distances. We would like to emphasize that the func- 
tion ilj{t) in the theorem is independent of the set Q and that there are no 
regularity assumptions on this set. If the set is regular, then there are few 
points at small distance from the boundary and the approximation is bad 
only on a small set. If the set is fractal, then there are many points at small 
distance from the boundary and the approximation is bad on a large set. 
In particular, this approximation is related to the Minkowski content of the 
boundary. See [S]. 

In the second section these approximations are applied to Erdos Turan 
estimates of irregularities of point distribution on a torus. In particular, 
following [16j and [19], we obtain explicit estimates for the discrepancy of 
sequences in lattices and in arithmetic progression, which improve and extend 
some results already in the literature. Here are some examples of the results 
in the second section (see, respectively. Corollaries [HI dSl and below, 
where these results are proved). 

Theorem 2 IfO<a<l,fiis the Lebesgue measure, Q is a measurable 
set m the torus T'^ = M°'/Z°', let 

M (a, n) = sup ({dist (x, dVl) < t}) 

t>o 

(this quantity is related to the Minkowski content of dVt). Also, if m^^'^ is 
a positive integer, let L{m) = m~^/'^7/' be the lattice points m~^^^ {gi, ...,gd) 



3 



with Qj = 0, 1, 
such m > 1, 



ra 



i/rf _ X ^ Then there exists a constant c such that for every 



sup 

Af(a,f7)<7 



m 



xGL{in) 



Theorem 3 Given e > 0, for almost every x in T'^ there exists a constant 
c such that for every m > 1, 



sup 



m) 



Theorem 4 Let X be a finite collection of hyperspaces o/R"' and let P{X) 
be the collection of all convex polyhedra in the torus T"' with diameter smaller 
than 1 and facets parallel to elements of X . Then there exists a constant c 
such that, given a prime number m, there exists a lattice point g = {gi, ■■■,gd) 
in if- with 1 < gj < m — 1, such that 



sup 

neP{X) 



IJ,{n) - m ^ V'xn {jm ^g) 



< cm ^ log'^l 



m] 



In particular, up to a logarithmic transgression, the discrepancy of a 
random arithmetic progression is comparable to the one of the lattice 

L(m) with the same number of points. On the other hand, while the first 
two theorems applied to polyhedra with a = 1 give the bound m~^l^, the 
third theorem gives the better bound m~^. This improves and extends a two 
dimensional result in [U Theorem 4D] , where it is proved that the discrepancy 
of m points with respect to polygons is dominated by log^'^^ (m). 

In the third section the results obtained on the Euclidean settings are ex- 
tended to compact Riemannian manifolds. We first consider approximations 
by linear combinations of eigenfunctions of the Laplace Beltrami operator on 
a compact Riemannian manifold. Then we state an analog of Erdos Turan 
estimates for irregularities of point distribution on manifolds and, inspired 
by [E], when the manifold is a compact Lie group or homogeneous space 
we consider point distributions generated by the action of a free group. In 
particular, using the Ramanujan bounds for eigenvalues of Hecke operators 
obtained in the above quoted paper, we prove the following result (see Corol- 
lary El below). 
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Theorem 5 IfAi = S0{3)/ S0{2) is the two dimensional sphere, ifH is the 
free group generated by rotations of angles arccos(— 3/5) around orthogonal 
axes, then there exists a constant c such that, if k is an integer and {(Tj}^^^ 
is an ordering of the elements in "H with length at most k, then for every x. 



This resuh has been proved in [T7] in the case of spherical caps, with a 
proof that relies on explicit estimates of Fourier coefficients. On the contrary, 
our result applies to more general domains. 

The authors wish to thank Luca Brandolini for several useful discussions 
on the subject of the paper. 

1 Approximation by entire functions 

The main result in this section is the following. 

Theorem 6 There exists a positive function ip{t) on t > with fast decay 
at infinity, ip(t) < c{a) (1 + for every a, such that for every measurable 
set fl in the Euclidean space M.'^ and R > 0, there exist entire functions A{x) 
and B{x) of exponential type R, satisfying 

A{x) < xn{x) < B{x), \B{x) - A{x)\ ^ i; {Rdist {x,dn)) . 

Proof. By a scaling argument, the statement for the set Q at the point 
X with functions of exponential type R, is equivalent to the statement for 
the set Rfl at the point Rx with functions of exponential type 1. Hence, 
it suffices to prove the theorem when R = 1. Let m(^) be a smooth radial 
function on R"^ with m(0 = if |^| > 1/2 and J^^mi^fd^ = 1. Then the 
convolution m * m{^) is a smooth radial function with m * m (0) = 1 and 
m(0 = if 1^1 > 1. Define 



This cumbersome definition will be clarified in a series of steps. 
Claim. The kernel K{x) is an entire function of finite exponential type, it is 
positive with mean 1 on R*^ and all its derivatives have fast decay at infinity, 
K{x) / dx^\ < c(l + |a;|)~" for every a and {3. 

All of this follows from the corresponding properties of the Fourier trans- 
form K{£^) = (l + 1^1^) ( + )/ ^ ^ m(^). Since this Fourier transform is 



m 



-1 




m* m{^) exp (27ri^ ■ x) d^. 
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smooth with compact support and it is 1 at the origin, the kernel has mean 
1 and all its derivatives have fast decay at infinity. Since the kernel is the 

convolution of the Fourier transform of (l + |^|^) ^ ''^ , which is up to a 
constant exp (— 27r |a;|), and the square of the Fourier transform of m (^), the 
kernel is positive. Finally, since K{^) = if |.^| > 1, by the Paley Wiener 
theorem the kernel is an entire function of finite exponential type, 

K{x + iy)^ / (l + ICO ^"^^^^^^ m*m{^)ex.p{2ni^- {x + iy))d^, 
\K{x + ly) I < (^j^^ (1 + lei') \m * m {i) \ d^^ exp (27r \y\) . 



Claim. Let 



I{t) = I K{x)dx. 



'{\x\>t} 

Then, for every t > 0, I{t + 1) > exp (-27r) I{t). 

Since K{x) is the convolution of cexp (— 27r \x\) and the positive function 
|m(a;)|^, one has 

K{x + y)= / cexp {—27: \x + y — z\) \fh{z)f dz 
> exp (— 27r ||/|) / cexp (—27r |a; — |m(2;) 1^ (i^; = exp (—27r i^r(a;). 



Hence 



I{t + 1)= / K{pd)p'^-^dpd'd 

jm=i}Jt+i 



r r+oo 

= / / K((T + l)^))(T + l)'^~Urd^ 
J{\^\=i}Jt 

r r+oo 

>exp(-27r) / / K {T^))T'^-^dTd^^exp{-2n)I{t). 



Claim. Define 

G{x) =/(dist {x,dn)). 

Then, for every x, 

Ixq {x) -K*xn{x)\<G{x). 
Since K{y) is positive with mean 1, 



\xn (x) -K*xn {x)\ = 



/ K{y){xn 



x) -Xn{x- y)) dy 



< [ K{y)dy = G{x). 
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Claim. Define 

F{x)^I (dist {x, dQ) /2) . 

Then, for every x, 

K*G{x)< 2F{x). 
Since dist {x — y, dO) > dist (x, dQ) — \y\, it follows that 



{\z\ > dist {x - y, dn)} Q {\z\ > dist {x, dn) - \y\} . 



Hence 



K*G{x)= K{y)K (z) dzdy 

JtR'^J {\z\>dist{x-y,dQ)} 



K {y) K {z) dzdy 

{|2/|<dist(a:,an)/2}^{|z|>dist(a;,aQ)/2} 



^/ / 

^•nj/|<dist(a:,an)/2}^{| 

+ / / K{y)K{z)dzdy 

J{\y\>dist{x,dn)/2}JM.d 

<2[ K(y)dy. 



Claim. Let 

= exp (-27r) / K (y) dy. 

Ah<i} 

Then for every x, 



K*G{x) >7-^G(x). 
Since dist {x — y, dQ) < dist (x, dfl) + \y\, it follows that 

{|^| > dist {x - y, dfl)} D {\z\ > dist (x, dQ) + \y\} . 



Hence 




K*G{x)= I K{y)K {z) dzdy 

' {\z\>dis\,{x-y,dn)} 



K iy) K {z) dzdy 

' {\y\<l}J {\z\>dis\,{x,dQ.)+l} 

> -f-^i (dist (x, an)) . 

To conclude the proof of the theorem, define H[x) = ■jG (x) and 

A{x)^K*xn {x)-K*H{x), 
B{x)^K*xn {x) + K*H{x). 



Since the kernel K(x) is an entire function of finite exponential type, then 
also the convolutions with this kernel are entire functions of finite exponential 
type. In particular, both A{x) and B{x) arc entire functions of exponential 
type not larger than the one of K{x). Moreover, by the above claims, 

Xa (x) - A{x) ^K*H{x)-{K*xn (x) - xn (x)) > 0, 
B{x) -xn{x)^K*H (x) - ixn (x) - K * xn (x)) > 0. 

Finally, 

B{x) -A{x)^2K*H (x) < 47F (x) = 47/ (dist (x, dQ) /2) . 

Hence, the theorem follows with ipij:) = 47/(^/2). ■ 

It follows from the proof of the theorem that for periodic sets with re- 
spect to the integer lattice Z*^, the above approximating entire functions are 
periodic too, hence they are trigonometric polynomials. 

CoroUciry 7 There exists a positive function ip{t) with fast decay at infinity, 
such that for every measurable set Q in the torus T'^ — R'^/Z'^ and R — 
0, 1, 2, there exist trigonometric polynomials A{x) and B{x) of degree R 
with 

A{x) < xq{x) < B{x), \B{x) -A{x)\ ^ {Rdist {x,dQ)) . 

Proof. The corollary follows immediately from the theorem. However, in 
order to clarify what follows, we write down explicitly the Fourier expansions 
of the trigonometric approximations. Let Q = Q + Z*^ be a Z'^-periodic set in 
R'^. As in the proof of the theorem, for every R > 0, let 

Kr (x) = ^'^^ {x + k))=J2^ i^/R) (27riA; • x) , 
//K(x)=exp(27r)f / K{y)dy\ [ K{y)dy. 
Then, 

A{x),B [x) = Kr {y) (xn {x-y)^ Hr {x - y)) dy 
= J2k {k/R) {xn {k) T Hr {k)) exp {27iik ■ x) . 
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Remark 8 As we have said, in the above theorem the approximation \B{x) — A{x)\ 
is controlled by 1 at points with distances 1/R from the boundary ofQ, while 
\B{x) — A{x)\ is essentially at larger distances. It follows from the in- 
equality of S.Bernstein between the maxima of an entire function and its 
derivatives that this approximation is essentially optimal. Indeed, if C{z) is 
an entire function of exponential type 1, then 

\C{x) - < sup \VC{z)\ \x-y\< 271 sup \C{z)\ \x - y\ . 

Hence, if A{x) < Xni^) < B{x) are entire functions of exponential type 1, if 
x is in VL and y is outside VL with \x — y\ < {Aiisn^i^^^'^d {\A{z) \ ,\B{z)\])~^ , 
then B{x) — A{y) > 1 and 

B{x) - A{x) = {B{x) - A{y)) - {A{x) - A{y)) > 1 - 1/2, 
B{y) - A{y) = {B{x) - A{y)) - {B{x) - B{y)) > 1 - 1/2. 



2 Discrepancy on the torus 



As advertised in the introduction, the approximation resuhs in the previous 
section have simple and straightforward apphcations to multi dimensional 
versions of the classical Erdos Turan inequality for discrepancy of point dis- 
tribution. In the sequel, for simplicity, the sets VL in the torus T'^ = M'^/Z'^ 
will be periodic sets in of the form = fi* + Z"', with Vt* in W'- with the 
property that dist {VL* ,VL* + k) > e > Q for every k e Z"^ - {0}. With this 
identification, the distance of a point x from dfl in T'^ is the distance of x 
from dn* + Z"^ in M*^. 

Theorem 9 If {^j}"=i o sequence of points in the torus, if Q is a mea- 
surable set with measure ^{fl), and if Hr{x) = 4"^?/' (2i?dist (x, with 
R > and ip (t) with fast decay at infinity, as in the proof of Corollary 
then 



i=i 



< 



^^(0)|+ E {\Xnm + \HR{k) 



0<\k\<R 



m exp {2Txik 



Proof. If A{x) < Xn{x) < B{x) are defined as in the proof of Corollary [71 
then 



A{x), B{x) =J2k (k/R) (xn (k) t Hr (k)) exp {27rik 



XI 
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Since 
then 



K{0 <l,K (0) = 1, K{0 = if lel > 1, and since xn (0) = /i (1^), 



/X (^]) - 5^ (xn (k) - Hr (A;)) m"!^ exp (27r2A; 



< 



0<|fc|<iJ 



m 

m~^y~^ exp (27riA; 



Similarly 



< 



0<|fc|<iJ 



m ""^^^ exp (27riA; ■ Xj) 



In order to apply the above theorem one has to estimate the exponential 
sums of point distributions and the Fourier transforms of domains. Motivated 
by the above result and by the definition of Fourier dimension of a measurable 
set, see e.g. Chapter 4.4 of [8], it is possible to introduce the classes of 
sets whose Fourier transform have a prescribed decay at infinity. Given a 
measurable set fl in the torus T'^, assume that for some < a < {d + l)/2 
there exists a constant c such that for every k — {0}, 



fd 



Xn (x) exp {—2Tiik ■ x) dx 



<c\k\-". 



Also, assume that for some < /3 < 1 and for every i? > 0, 

c\k\-'^ if < \k\ < R, 



4^ (-Rdist (x, dVL)) exp {—2nik ■ x) dx 



< 



cR-P if = 0. 



Denote by F (a, /3, Vt) the smallest constant c for which the above two 
inequalities hold. It turns out that the Fourier transform of ip {R dist (x, dQ)) 
is in some sense dominated by the one of xn{x)- In particular, in a number 
of cases the second inequality is a consequence of the first. Anyhow, in many 
cases it is possible to give quite precise estimates for the constants F (a, (3, fl). 
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Remark 10 If the first of the above inequalities holds for some a > {d+l)/2, 
or if the second inequality holds for some P > 1, then either D, orT'^ — Q has 
measure zero. It suffices to show this when {d + l)/2 < a < d/2 + 1. For 
every y, 

/ \Xii{.x + y) -xn{.x)\^dx = ^ \{exp{27rik ■ y) - 1) xn{k)f 

V 0<lfe|<|y|-i \k\>\y\-' 

<cF(a,/3,0)|t/r-^ 

Hence for every y and every n> \y\, 

/ \xn{x + y) - xUx)\ dx 

n „ 

< / \xn{x + jn~^y) -xn{x + {j - i)n~^y) \ dx 

n „ 

= X] / + j^~^y) - (a; + {j - ^ dx 

j \xq {x + n~^y) - xn {x)\'^ dx < cF{a,l3,Q) \yf°''''^n' 



2a-d d+l-2a 



This converges to zero when n diverges. Hence, for every translation Q — y = 
Q up to sets with measure zero. Similarly, if 13 > 1, then either Q orT'^ — Q 
has measure zero. To see this, it suffices to estimate the modulus of continuity 
of the function 

G{y)= \xaix + y)-xnix)\dx. 

Indeed, 

\Giy)-G{z)\ 

< / \xn{x + y) - Xn{x + z)\dx < ii{{dist{x,dn) <\y - z\}) 
< [ inf {V^ {t)}] I il}{\y- z\-^ dist (x, dO)) dx < cF {a, /3, 0) \y ■ 

\0<t<l J Jjd 

Finally, if (3 > 1, then G{y) is identically zero. 



1/3 
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Remark 11 The Fourier coefficients of ip {R dist {x,dfl)) on the torus T'^ 
can be evaluated by an integration over any set that tiles M'^ via TL'^ . In 
particular, one can integrate on the set of points Q in M°' for which the 
distance from dQ + Z*^ is realized precisely by the distance from dQ. Since 
I V dist {x,dQ)\ = 1 when dist {x,dQ) ^ 0, if dQ has measure zero, then the 
coarea formula gives 



ip {R dist (x, dQ)) exp {—2ni^ ■ x) dx 



Qn{dist(x,an)<e} 
eR 



pen 
^0 



exp [—2711^ ■ x) dx dt. 



' {dist{x,dn)=R-'^t} 

Moreover, since Q has measure 1, 



ip (-Rdist (x, dQ)) exp (— 27ri^ ■ x) dx 



Qn{dist(2:,ar2)>e} 



< sup \ij{t)\ 

t>eR 



In fact, one can eliminate the term sup^^^^ {ip {t)\ by integrating over all 
Q rather than Q (1 {dist {x,dQ) < e}. However, in order to keep control of 
the level sets {dist (x, dQ) = t} it may be convenient to restrict to t small. 

Remark 12 It follows from classical estimates on oscillatory integrals with 
non degenerate critical points that if a convex body has smooth boundary with 
positive Gauss curvature, then 



/ exp {—2Tci^ ■ x) dx 
Jan 

exp {—2'KiC, ■ x) dx 



<d-l)/2 



^(d+l)/2 



The constant c can be bounded in terms of the smoothness and the min- 
imum of the curvature of the boundary. See e.g. /TOj [77|. 17^ \2^. From 
this estimates, with {dist (x, dQ) = t} instead of dQ, it easily follows that 
F {{d + l)/2, 1, Q) is finite. If the curvature vanishes of some order at some 
point, then the Fourier transform in directions normal to these points has a 
worse decay at infinity. If some part of the boundary is completely flat, then 
one can guarantee only a decay of order one. 

Remark 13 If < a < I, define 



M (a, Q) = sup ({dist (x, dQ) < t}) . 
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This is related to the upper Minkowski content of dQ, defined by 
lim sup ({dist {x, dil) < t}) . 

i-i-0+ 

However, these two quantities can be quite different. In particular, if dQ 
contains a point p, then M {a,Q) > fi {{dist {x , p) < 1}). If M {a,Q) is 
finite, then dQ has Minkowski dimension at most d — a. The domain is 
regular if a = 1 and it is fractal if < a < 1. The definition makes sense 
also when a > 1, but in this case either Q orT'^ — Q has measure zero. The 
proof is as in Remark [73 It is well known that the decay of the Fourier 
transform of a domain can be controlled in terms of these quantities. Indeed, 

Xn{x) exp{—2'Kik ■ x)dx 

exp(— 27riA; • [x — 2^^ \k\^'^ k^)dx 

= / (xnix) -xn{x + 2'^ Ikf^ k)) exp{-27rtk ■ x)dx. 

Then, 

Xn{,x) exp(— 27rzfc ■ x)dx 

<2-^ I \xn{x + 2-^\k\-^k)-xn{x)\dx 
< 2-V {dist (x, dn) < 2-1 \k\-^] < 2-"-iM (a, fi) l^"" . 
Moreover, since ip [t) is positive and has fast decay at infinity, 

ip (-Rdist (x, dVt)) exp {—2TTik ■ x) dx 
< [ ip {R dist {x,dn))dx<cM {a, n)R-". 

In particular, since i?"" < in the range < \k\ < R, it follows that 

F{a,a,n) < cM (a,fi). 

In the following, the above theorem and remarks will be applied to the 
study of discrepancies of lattices m~^^^'Z'^ and arithmetic progressions {jx} Jl^ 
in the torus M'^/Z'^, with multiples modulo Z'^. In particular, we will prove 
Theorems [21 [3] and [H stated in the Introduction. 
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Corollary 14 There exists a constant c with the following property. Let 
m}^'^ he a positive integer and let L{m) = m~^/°'Z°' be the lattice of points 
m-i/°' (gi, ...,gd) with < gj < m}!"^ - 1. Then 



x(iL{m) 



< cF{a,/3,n)m' 



-/3/(d+/3-a) 



Theorem [2] follows from the above Corollary by setting /3 = a, replacing 
F{a, a, Q) with M(a, Q) and observing that the constant c does not depend 
on the set fl. 

Proof. The sum of a geometric progressions gives 



m 



exp {27Tik ■ x) 



x(^L{m) 



n=i y j=0 J 



if A; ^ m^/'^Z'^, 
if A; e m^/^'Z'^. 



Hence, by Theorem [9] and the definition of F{a, 

<F{a,/3,n) R~^ + 2 

y 0<|ml/dfc|<H 

Then the choice R = rn}^^^~^^~°'^ gives the desired estimate. ■ 



Xn (x) 

x(^L{m) 



\ 



m 



Corollary 15 Given e > 0, for almost every x in T'^ there exists a constant 
c such that for every m > 1 and every fl, 



As before, Theorem [3] follows from the above Corollary by setting /3 = a, 
replacing F{a, a, Q) with M{a, Q) and observing that the constant c does 
not depend on the set Q. 

Proof. Denoting by the distance of t to the nearest integer. 



m ""^^^ exp (2TTijx ■ k) 



sin (Tcmk ■ x) 



msin (tt/c ■ x) 



< min {1,1/ (2m ||A; ■ . 
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Hence, by Theorem [H] and the definition of F{a, 



< j i?-^ + 2 \k\-"mm{l,2-^m-^\\k-xr^} 

0<\k\<R 
y 0<|A:|<_R 

Finally, in [12] it is proved that for almost every x there exist a c which 
depends on x, such that for every R > 0, 

J2 |A;r'^||A;-x|r' <clog'^+'+"(l + i?). 

0<|fc|<i?, 

The desired result follows choosing 

■ 

In particular, the discrepancy of a random sequence {jx}^^^ with respect 
to any domain is dominated by cF {(d + 1) / 2,1, n)m~'^^^^^^Hog^~^^{m) 
and, by Remark [121 when Q is convex with smooth boundary and positive 
Gauss curvature, then F {{d + l)/2, 1, fl) is finite. Similarly, by Remark [T5| 
the discrepancy is dominated by cF (1, 1, Q) m~^^'^ \Qg{<i+^+£)/<^ [rn), and when 
the boundary of the domain is (i — 1 dimensional, then F (1, 1, i7) is finite. 
These results should be compared with well known upper and lower estimates 
of the discrepancy with respect to convex regions due to W.M. Schmidt and 
J. Beck. See e.g. 0, Theorem 15 and Corollaries 17B, 18C and 19F]. In par- 
ticular, for any given compact convex set there exists an infinite sequence 
such that the discrepancy is bounded above by cm"'-'^"'"^^/^^'^^ log^''^ The 
definition of Fourier dimension does not capture polyhedra, except in the 
trivial case a < 1. Indeed the decay of Fourier transforms of polyhedra 
is not isotropic nor homogeneous. Anyhow, these Fourier transforms can 
be computed explicitly and estimated quite precisely and, using these esti- 
mates, one can give bounds for the discrepancy which are better than the 
ones obtained above. 

Lemma 16 If Q is a polyhedron in R'^ with diameter X, then, 

P d 

exp{-2TTi^-x)dx <2 ^ JJmin |a, (27r |Pf^(j)^|)"^| . 
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The sum is taken over all possible decreasing chains of j dimensional faces 
Q(j) of D, and Pny) is the orthogonal projection on the j dimensional sub- 
space parallel to fl {j). 

Proof. The Fourier transform of a polyhedron can be computed exphcitly, 
but here we are only interested in precise estimates of its size with control on 
all constants involved. If Q (j) is a j dimensional face of Q with Pn{j)C 
and if n (x) is the outgoing normal to the boundary dfl (j) at the point x, 
the divergence theorem gives 

exp (— 27ri(^ ■ x) dx = f - — ^— ^ [Hi^ ^^^p (^_2^j^ . x'j dx. 

Moreover, if A is the diameter and ii (^^(i)) is the j dimensional measure of 
fl{j), one always has the trivial estimate 



exp (— 27ri^ • x) dx 



</i(Q(j))<A^'. 



Hence, if Q (j — 1) are the j — 1 dimensional faces of (j). 



exp (— 27ri^ • x) dx 



< 



E 



27r \P^i(j)^\ no-i)cno) 



/ 



exp (— 27ri^ • x) dx 



if 27r|Pf,(,)e| < 1/A, 
if27r|Pf,(,)e|>l/A. 



Iterating, one can decompose the integral over Q = fl[d) into a sum of 
integrals over chains of faces n{d) D n{d — 1) D ... and this gives 



exp (— 27ri^ • x) dx 



n{d)DQ{d-l)D...DQ{s) j=s+l 
d 

+2 Yi n(2^i^^o)^i)"'- 

n{d)Dn{d-i)D...Dn(i) j=i 



The first sum is over all chains of faces with (27r < A for 1 < s < 

j < d and (27r |i-*n(s)C|) ^ > while the second sum is over all chains with 



16 



(271 \ Pn{j)^\) ^ < A for all 1 < j < d. Finally, in the first sum, 

d d 
j=s+l j=l 

-1 



Indeed, since |-PQ(j)C| is increasing in j, the terms min|A, (27r 

are equal to A when 1 < j < s and equal to (27r |-Pn(j)^|) ^ when s < j < d. 
Similarly, in the second sum, 

d d 
3=1 3=1 

Observe that when ^ = the formula gives A'^, while the exact value of the 
integral is the volume of Q. ■ 

Lemma 17 Let Q be a convex polyhedron in M'* with diameter A. For any 
j = 1,2, ...,d — 1, let {A (j)} be the collection of all j dimensional subspaces 
which are intersections of a number of subspaces parallel to the faces of Q. 
Finally, let 'ip{t) be a function on < t < +oo with fast decay at infinity. 
Then there exists a positive constant c, which depends on d and'ip{t), but not 
on fl, such that for every R> 0, 



/ 



-0 (i?dist {x, dfl)) exp {—2ni^ ■ x) dx 
E ^'-'riniin{A,(27r|P^(,)e|)~'}. 



d-l 



3=0A{j)^...DA{l) k=l 

When j = the inner sum of products is intended to be the number of 
vertices of the polyhedron, when 1 < j < d — 1 the inner sum is taken over 
all possible decreasing chains of j dimensional subspaces {A{j)} and Pa^j) is 
the orthogonal projection on A (j) . 

Proof. Since |Vdist {x, d^) \ = 1, the coarea formula gives 

ijj (it!dist {x, do,)) exp (— 27ri^ • x) dx 



+ CXJ 



exp (— 27ri^ ■ x) dx\ i/j (Rt) dt. 

{dist{x,dn)=t} / 

We consider separately the level sets inside and outside ^2. The level 
sets {dist {x, dfl) = t} nfl are polyhedra with diameter at most A and faces 
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parallel to some {A{j)}. We emphasize that some of the {A{j)} may not 
be parallel to any {Q (j)}- Anyhow, as in the previous lemma, 







exp {—2ni^ ■ x) dx\ i/j (Rt) dt 



{dist(a;,9f2)=t}nO 
oo 



^ E n \^ W\ dt) ilmin {a, {2n \PAij)^\y'} ■ 

A{d-1)D...DA{1) V -^0 / j=l 

As in Steiner formula, the outer level sets {dist {x, dfl) = t} — fl are union 
of sums of j dimensional faces {Q (j)} and portions of — j — 1 dimensional 
spherical surfaces of radius t. Hence, if uj{d — j — 1) denotes the measure of 
the d — j — 1 dimensional spherical surface of unit radius. 



+00 



exp (— 27ri^ ■ x)dx] ijj (Rt) dt 



{dist{x,dn)=t}-n 



< 



j-i; 



+00 



\iP{Rt)\t''-^-''dt 1 X 



j2Ud 

j=0 ^ 

j 

Cl{j)D...Dn{l)k=l 

Observe that when = the formula gives c (^R~^ + X^~^R~^^, with a con- 
stant c independent of the polyhedron. ■ 

Theorem 18 Given a finite collection of d — 1 dimensional hyperspaces X 
inM.'^, letP{X) he the collection of all convex polyhedra with diameter smaller 
than 1 — e and facets parallel to elements of X . If {A (d) D ... D A (1)} is 
the collection of all possible decreasing chains of j dimensional subspaces 
obtained by intersection ofW^ and a number of hyperplanes in X , define 



*(0 



a 

J2 llmm[l,{27r\PAU)^\r'} 



A(d)D...DA(l) j=l 



Finally, if {xj}"^^-^^ is a sequence of points in the torus, define 



m exp (27ri^ ■ xj) 



Then, there exists a positive constant c, which depends only on the space 
dimension d and the upper bound 1 — e of diameters of polyhedra, such that 
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for every R > 0, 



sup 



0<\k\<R 



Proof. This follows from Theorem [9] and the previous lemmas. It suffices to 
replace X < 1 — e with 1 and R^~'^ with in the range < \k\ < R. ■ 

For example, if X is the collection of hyperplanes {xj = 0}, then P{X) 
is the collection of all boxes {aj < Xj < bj} with diameter smaller than 1 — e 
and the above is an estimate of discrepancy with respect to boxes. 

Corollary 19 Given a finite collection of hyperspaces X and a prime num- 
ber m, there exists a lattice point g = {gi, ■■■gd) in with 1 < gj < m — 1, 
such that 



sup 



- m ^^Xn {jm ^g) 



< cm ^ log'^(m) 



The constant c depends only on the dimension and on the cardinality of the 
set of chains of subspaces A (d) D . . . D A (1) generated by X . 

Theorem H] follows from the above Corollary, along with the observation 
that a polyhedron in the torus with facets parallel to elements of X and 
diameter smaller than 1 can be seen as the union of a finite number (which 
depends only on X) of polyhedra with facets parallel to elements of X and 
diameter smaller than, say, 1/2. 
Proof. A preliminary result is needed. 

Claim. There is a constant c, depending only on the dimension d, such that, 
for any decreasing chain of subspaces A (d) D . . . D A (1) and for any R > 0, 

d 

J2 Umin {^1, {2n \PAU)k\y'] < clog'' {2 + R). 

l<\k\<Rj=l 

Indeed, let {ai, . . . , a^} be an orthonormal basis of M'^ such that {ai, . . . , Oj} 
generates A{j) for j = 1, d. Since |PA(j)fc| > \k ■ aj\, it suffices to estimate 

d 

l<\k\<Rj=l 

Of course the idea is to replace the sum over a discrete variable with an 
integral over a continuous variable. Let Q{k) be the cube centered at the 
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integer point k with sides of length 1 parallel to the orthogonal axes {aj} and 
let P{k) be the cube centered at k with sides 3. Since \x ■ aj\ is the distance 
of X from the hyperplane of equation x ■ Oj = 0, if this hyperplane crosses 



P{k), then \x ■ aj\ < 2\fd for any x in Q{k) and therefore 

mm = mm \x ■ aj\ > . 



If this hyperplane does not cross P{k), then for any x in Q{k), 



mm 





X ■ 






\k- 






X 


aj\ 



< min < 1, 



\k ■ Qjl 


+ Vd/2 


1 /l ' Qjj 





> < min |l, ^1 + Vd/2) I X ■ ail > . 



In the overall, for any k and any x G Q{k), 

d d 

JJ^min (l, \k ■ a^P^) < JJ^min ( l,2v^| X ■ aA 



Hence, 



J2 Ylmm{l,\k-aj\ ^} 

l<\k\<Rj=l 
d 

dx 



< ^ / JJ min 1 1, 2 Ix ■ a^ r^ I 

d 

TTmin|l,2v^ 

V2},=1 ^ 



< / II min X ■ aA > dx 



{\x-aj |<-R+1 



2 



rR+l/2 ^ . Y 

J min|l,2Vrfr^|dt j < c log"* (2 + i?) . 



Now comes the proof of the Corollary. The sum of a geometric progression 
gives 

E"^ /r. ■ ■ -^ 7\ smiiTg-k) , ., _i ,n 

exp [27rtjm g ■ k) = — exp [mim + l)m g ■ k) . 
^ ' msm[Tim~'^g ■ k) ^ ' 

This exponential sum is or 1 according to ■ 7^ (modm) or g ■ k = 
(modm). Hence, by Theorem [TS] with R = m, 



sup 

neP{x) 



- m ^^Xn {jm ^g) 



< c j m'^ + ^ $ (A;) 

0<|fc|<m, g-k=0 (modm) 
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The heuristic behind the existence of good lattice points with the de- 
sired properties is that the ratio of fc's which satisfy the congruence g ■ k = 
(modm) is and the sum over the fc's with g ■ k = (modm) is m~^ 
times the sum over all k. This heuristic principle can be made rigorous by 
an averaging procedure, as in Theorem 5.7 of |l6j. In order to satisfy the 
congruence giki + ... + g^kd = (modm), if /cj 7^ (modm), then one can 
take gj arbitrary for j 7^ i, the remaining gi being uniquely determined in 
the residue class. Hence, by the above claim, 

l<gj<m— 1 \ 0<|fc|<m, g-fc=0 (modm) 

= E '^(^) (("^-1)"' E 1 

0<|fc|<m \ g-k=0 (modm) 

0<\k\<m 

Observe that the constant c is the product of the constant in the above 
claim and the cardinality of the set of chains of subspaces generated by X. 
In particular, there exists g such that 



E 



$ (k) < cm log (m). 

0<|A:|<m, g-k=0 (modm) 



Corollary 20 Given a finite collection of hyperspaces X and e > 0, then for 
every e > and almost every x in T'^ there exists a positive constant c such 
that for every m > 1, 



sup 

neP{x) 



< cm-Mog'^+^+'l 



m 



Proof. As in Corollary [151 this follows from Theorem [18] and an adaptation 
of[l9]. ■ 

The above corollaries improve and extend a two dimensional result in 
[H Theorem 4D], where it is proved that the discrepancy of m points with 

m). They should also be 



4+£ 



respect to polygons is dominated by m ^ log' 
compared with the discrepancy of lattice points in [3] . 
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3 Approximation and discrepancy on mani- 
folds 



Let Aihe a smooth d dimensional compact manifold without boundary, with 
Riemannian distance dist(a:, y) and measure /i normalized so that /i (Ai) = 1. 
The Laplace Beltrami operator A on has eigenvalues {A^}, counted with 
appropriate multiplicity, and a complete orthonormal system of eigenfunc- 
tions {(p\{x)}. To every function in L'^{J^,dfi) one can associate a Fourier 
transform and a Fourier series. 



Fourier series on compact Lie groups and symmetric spaces are examples. 
In particular, the eigenfunctions of the Laplace operator on the torus are 
trigonometric functions and eigenfunction expansions are classical Fourier 
series. Similarly, eigenfunctions of the Laplace operator on the surface of a 
sphere are homogeneous harmonic polynomials and eigenfunction expansions 
are spherical harmonic expansions. In the setting of manifolds, an analog of 
trigonometric polynomials is given by finite linear combinations of eigenfunc- 
tions ^^^cxipx^x). Indeed it can be shown that there is a close relation be- 
tween approximation by functions of exponential type and by eigenfunctions. 
See for example [6j . The following generalizes Theorem and Corollary [71 

Theorem 21 Given a > 0, there exists /3 > such that, for every do- 
main Q in Ai and R > 0, there exist linear combinations A{x) and B{x) of 
eigenfunctions with eigenvalues at most , satisfying 



It is likely that a slightly more precise result holds, with a rapid decay 
instead of a polynomial decay (1 + fZdist (a;, 9^2))"°. However, the exponent 
a can be arbitrarily large and this suffices for our applications. 
Proof. The proof of this theorem is similar to the proof of Theorem [H] and it 
is based on suitable approximations of the identity adapted to the manifold, 
analogous to the convolution kernels in the Euclidean spaces: 




A 



A{x) < xn{x) < B{x), \B{x)-A{x)\ ^ (3 (1 + R dist {x , dn)) 



—a 




\<R 



Kr{x, y)\<c (a) R\l + R dist (a;, y)) 



a 
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Moreover, these kernels are positive up to a negligible error. The con- 
struction of such kernels on Lie groups and symmetric spaces is well known 
and in these cases it is possible to obtain positivity. Indeed, if a kernel has 
good decay and finite spectrum, then also its square has good decay and 
finite spectrum and a suitable normalization has mean one. We do not know 
whether positivity can be achieved in our general setting, however in the 
sequel almost positivity will suffice. Given m (^) as in the proof of Theorem 
[6] and 

Mlel) = (l + l^l')"^'^'^^'rn*m(0, 

define 

KR{x,y) = (R'^X) (px{x)(pxiy). 

A 

It is possible to prove that this kernel has an asymptotic expansion with 
Euclidean main term /2'^A'(i?dist(a;, y)) and suitable control on the remain- 
der. Although the details are not completely trivial, the techniques can be 
found in Chapter XII of [2T], or in [3]. Finally, define 



A{x) 
B{x) 



Hr{x) = /3(l + i?dist {x,dn)y, 

Kr{x, y) ixniy) - Haiy)) dfi{y), 



M 



Kr{x, y) {xn{y) + HR{y)) d^{y). 



M 



Then, as in the proof of Theorem [HI it is possible to show that for some 
/9 > independent of VL and R these functions satisfy the required properties. 



Theorem 22 For every sequence of points {xj}^^^ 
and R > 0, if Hr{x) is defined as in the proof of Theorem 



and domain Q in A4 
then 



< 



Hr{0)\+ J2 {\XnW\ + \HR{X) 



0<X<R 



m 

mT^^^Lfx {xj 



Proof. This proof is completely analogous to the one of Theorem [91 ■ 

Of course, the interest of the above result arises when one is able to 
exhibit point distributions with rn~^J2f=i^\{^j) suitably small. Inspired 
by [IZ] on the problem of distributing points on a sphere, we now consider 
point distributions generated by the action of a free group on a homogeneous 
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space. Let Q he a compact Lie group, /C a closed subgroup, Ai = Q/JC a 
homogeneous space of dimension d with normahzed invariant measure /x. Let 
"H be a finitely generated free subgroup in Q and assume that the action of 
"H on is free. Given a positive integer k, let {cjjjlj^ be an ordering of the 
elements in "H with length at most k. For every function f{x) on Ai, define 



Tf{x) 



m 



This operator is self adjoint with norm 1, hence all its eigenvalues have 
modulus at most 1. Indeed, 1 is an eigenvalue and the constants are eigen- 
functions. In the following, we shall be interested in cases where all other 
non constant eigenfunctions have eigenvalues much smaller than 1. For this 
reason, define p(m) as the supremum of the eigenvalues with non constant 
eigenfunctions, 

p(m) = sup |z/| . 

Moreover, as before, define M (5, Q) = sup^^o t~^fi ({dist (x, dQ) < t}). 



Theorem 23 There exists a positive constant c such that for every point x 
in M. and R> 0, 



sup 

Af(<5,n)<7 



Proof. Since the operators T and A commute, they have a common or- 
thonormal system of eigenfunctions, A(px{x) = X^ipxi^x) andTipx{x) = T{\)ipx{x). 
The assumption in the theorem is precisely that |T(A)| < p{m) if A 7^ 0. 
Hence, by Theorem [22| 



i=i 



< 



Hr{0) +p{m) J2 (\XnW\+ Hni\))\vx{p)\- 



0<X<R 
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Since Hr {x) = (1 + i?dist (x, dQ))'", then 

^r(0) < /3 / (1 + i?dist (x, c//i(a;) 
Jm 



+00 



< /3 I /i ({dist (x, dn) < R-^}) + ^2-"> ({dist (x, < 

j=0 

< /3 [ 1 + 2^^2(^-")M M (5, n) R-^. 



j=0 

Similarly, by Cauchy and Bessel inequalities, 



1/2 f >^ 1/2 

2 



</3 
One also gets 



J2 \HRix)\ i^A {p)\ < [yI I E 

</3|^(l + /2dist(x,ai]))-'"rf/i(x)| ^ |5^|</^a(p)| 

+OC >| 1/2 J- ^ 1/2 

j=0 J Vx<R ) 



1/2 



-<5/2 



E IXn(A)||^A(p)| 



0<A<_R, 

N 1/2 J. 1/2 

.0<A<1 j lo<A<l 

1/2 , . 1/2 



[log2(R)] f I f 

+ E E E 

A:=0 [^A>2* j U<2fc+i 



I^a(p)| 



2 



If A2k{x) < x^{^) ^ B2k{x) are the approximating functions in Theorem 
[2T] corresponding to i? = 2'^, then, 

ElXn(A)|'</ |xn(a;)-A.(x)|2d/i(^) 

A>2^ -^-^ 

</32/ (l + 2'=dist(x,91]))"'"ci/i(x) 
< /^^ ( 1 + 2'5^2(''-2a)n M (S, n) 2-^\ 



j=0 
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By classical bounds on the spectral function of an elliptic operator, see 
e.g. Theorem 17.5.3 in [T3] . 



\<R 

Hence, if a is large enough. 



< cM (5, fi) R-^ + c^M{5,n)R^'^-^^'^p{ 



m). 



Finally, observe that a/M {6, Q) < cM{6,Q), since M [6, Q) is bounded 
below as one sees putting t = 1 in the definition of this constant. ■ 

The following corollary is Theorem [5] in the Introduction, and it has been 
proved in |17J in the case of spherical caps. 

Corollary 24 If M. = SO{3)/SO{2) is the two dimensional sphere, ifH is 
the free group generated by rotations of angles arccos(— 3/5) around orthog- 
onal axes, then there exists a constant c such that, if k is an integer and 
{'^j}^=i ^■^ '^'^ ordering of the elements in % with length at most k, then for 
every x, 



X 



< cM {6, Q) m-^/(2+5) iog2V(2+5)(^)_ 



Proof. The eigenvalues of the operator T satisfy the Ramanujan bounds 
p(m) = sup l^l < cm~^^^ \og{m) . 

Hence, choosing R = m}/^'^^^^ log"^''*-'^"'''^^ (m) in the above theorem, 

inf {R-' + R^^-'^/^m-^/Hogim)] < cm-'/'-''+'hog''/^''+'\m). 
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